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Abstract. Could the discrepancies found in the determination of mass 
in clusters of galaxies, from gravitational lensing data and from X-rays 
observations, be consequence of the standard description of the ICM, in 
which it is assumed hydrostatic equilibrium maintained by thermal pres- 
sure? In analogy to the interstellar medium of the Galaxy, it is expected 
a non-thermal term of pressure, which contains contributions of magnetic 
fields. We follow the evolution of the ICM, considering a term of mag- 
netic pressure, aiming at answering the question whether or not these 
discrepancies can be explained via non-thermal terms of pressure. Our 
results suggest that the magnetic pressure could only affect the dynamics 
of the ICM on scales as small as <^ lkpc. These results are compared 
to the observations of large and small scale magnetic fields and we are 
successful at reproducing the available data. 



1. Introduction 

Since the work of Loeb & Mao (1994), the possibility of explaining the discrep- 
ancies on mass determinations, found by Miralda-Escude & Babul (1994), via 
non-thermal pressure support has been widely discussed. The discrepancy arises 
from the two most promising techniques to obtain clusters of galaxies masses. 
On one hand, the determination of masses in clusters of galaxies, via X-ray 
data, is based on the hypothesis that the ICM is in hydrostatic equilibrium with 
the gravitational potential, using the radial profiles of density and temperature 
(Nulsen & Bohringer 1995). On the other hand, gravitational lensing measures 
the projected surface density of matter, a method which makes no assumptions 
on the dynamical state of the gravitating matter (Miralda-Escude &: Babul 1994; 
Smail et al. 1997). 

In clusters with diffuse radio emission X-ray observations can give a lower 
limit to the strength of the magnetic field (the 3K background photons scattering 
off the relativistic electrons produces the diffuse X-ray emission). Typically, this 
limit is B > 0.1 /uG (Rephaeli et al. 1987) on scales of ~ 1 Mpc. Such a kind 
of detection of clusters magnetic fields leads, using ROSAT PSPC data and also 
327MHz radio map of Abell 85, to an estimate of (0.95 ± .10) /xG (Bagchi et al. 
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1998). In the case of Faraday rotation the information obtained is the upper limit 
on the intensity of the field, and the measured values are (RM < 100 rad/m 2 ), 
which is more or less consistent with a intracluster field of B ~ 1 p,G, with a 
coherence length of Ib < 10 kpc. This strength of the magnetic field corresponds 
to a ratio of magnetic to gas pressure ps/Pgas < 10~ 3 , implying that B does 
not influence the cluster dynamics (at least on large scales) . In inner regions the 
magnetic fields are expected to be amplified due to the gas compression (Soker & 
Sarazin 1990). For frozen-in fields and homogeneous and spherically symmetric 
inflow, B oc r~ l and RM oc r _1 , (pb <x r~ 2 whereas gas pressure increases). 
Very strong Faraday rotations were observed (RM ~ 4000 rad/m 2 ) implying 
B > 10 nG at l B ~ 1 kpc (Taylor & Perley 1993; Ge & Owen 1993, 1994). 

2. Evolution of the ICM with Magnetic Pressure 

Using a spherically symmetric finite-difference scheme Eulerian code, the evolu- 
tion of the intracluster gas is obtained by solving the hydrodynamic equations 
of mass, momentum and energy conservation (see Friaca 1993), coupled to the 
state equation for a fully ionized gas with 10% helium by number. The mass 
distribution, M(r), is due to the contribution of the X-rays emitting gas plus the 
cluster collisionless matter (which is the sum of the contributions of galaxies and 
dark matter - the latter being dominant) following p c i(r) = p c (l + r 2 /a 2 )~ 3 / 2 , p c 
and a (the cluster core radius) are related to the line-of-sight velocity dispersion, 
a, by 9cr 2 = AirGa 2 p c . The total pressure pt is the sum of thermal and magnetic 
pressure, e.g. pt = P + Pb- The constraints to the magnetic pressure come from 
observations, from which pb = B 2 /8it ~ 4 x 10~ 14 erg cm _3 s~l (cf. Bagchi et 
al. 1998) for a diffuse field located at ~ TOOfc^kpc from the cluster center. 

The initial conditions for the gas are an isothermal atmosphere (To = 10 7 K) 
with 30% solar abundance and density distribution following that of the cluster 
dark matter. The evolution is followed until the age of 14 Gyr. We assume: 
frozen-in field; spherical symmetry for the flow and the cluster itself; and that 
at r > r c (the cooling radius), the magnetic field is isotropic, i.e., B 2 = B 2 /2 = 
B 2 /3 and l r = l t = l (where B r and B t are the radial and transversal components 
of the magnetic field B and l r and l t are the coherence length of the large-scale 
field in the radial and transverse directions). In order to calculate B r and Bt for 
r < r c we modified the calculation of the magnetic field of Soker & Sarazin (1990) 
by considering an inhomogeneous cooling flow (i.e. Mj ^ M varies with r). 
Therefore, the two components of the field are then given by D / Dt(B 2 .r 4 M~ 2 ) = 
and D / Dt(B 2 r 2 u 2 M~ l ) = 0. In our models it is admitted that the reference 
radius is the cooling radius r c . In fact, we modify the geometry of the field 
when and where the cooling time comes to be less than 10 10 yr. Therefore, our 
condition to assume a non-isotropic field is t coo = 3kBT/2/j,mH^(T)p ^ 10 10 yr. 

3. Models and Results 

There are four parameters to consider in each one of the models: a = 1000 km/s, 
the cluster velocity dispersion; po= 1-5 x 10 _28 g cm" 3 , the initial average mass 
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Figure 1. Evolution of the density and temperature (left) and mag- 
netic field strength (right) profiles. Curves represent early and late 
stages of the ICM evolution, as labeled, for the model A. 

density of the gas; a = 250 kpc, the cluster core radius; and (3q = 10~ 2 (model 
A), 10 -3 (model B), the initial magnetic to thermal pressure ratio. 

First of all, the evolution we follow here is characteristic of cooling flow 
clusters and in this scenario we discuss the evolution of the basic thermody- 
namics parameters. Considering the overall characteristics of our models, we 
compare the present models with Peres et al. (1998) deprojection results (based 
on ROSAT observations), pointing out that the central cooling time here adopted 
as our cooling flow criterion, e.g. t coo ^ 10 10 yr, is typical for a fraction between 
70% and 90% of their sample. This allow us conclude that our models, which 
present cooling flows since the cluster has the age of ~ 7 — 9 Gyr, are typical for 
their sample. 

Figures 1 shows the evolution of density, temperature and magnetic field 
strength, for model A, from which the presence of the cooling flow on later stages 
of the ICM evolution and at inner regions, is remarkable, if one notices the steep 
gradients of these quantities. We chose two values of magnetic field strength 
for the comparison of the results, on small and large scales (see figure). Our 
results for the magnetic field strength and also for the pressure, on large and 
small scales, are in agreement to the observed ones. 

Obviously the magnetic pressure (Figure 2) is compatible with the magnetic 
field intensities and may be compared to the values determined by, for instance, 
Bagchi et al. (1998), p B = B 2 /8ir ~ 4 x 10" 14 erg cm^s" 1 , at scales of 700 kpc, 
in the present time. From the analysis of the magnetic pressures expected from 
our models it is clear that they agree, as well as the magnetic field strength, 
with the observations. 

The present models are in many aspects similar to the one of Soker & 
Sarazin (1990). However there are two important differences between our model 
and theirs: i) they take into account only small-scale magnetic field effects; and 
ii) they consider the magnetic field isotropic even in the inner regions of the 
cooling flow. As a matter of fact the magnetic pressure reaches equipartition 
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Figure 2. Evolution of the magnetic (dashed line) and thermal (full 
line) pressures profiles on late stages of the ICM evolution for the model 
A (/% = 1(T 2 , left) and model B (ft = 1(T 3 , right). 



only at radius as small as ^ 1 kpc (model A) or ;> 0.5 kpc (model B), because 
the central increase of the (3 ratio is moderate in our model. Our more realistic 
description of the field geometry is crucial. This implies that the effect of the 
magnetic pressure on the total pressure of the intracluster medium, even on 
regions as inner as few kpc, is small. Evolutive models for the intracluster 
medium, with a realistic calculation of the geometry and amplification of the 
magnetic fields, like the one presented here, indicate that magnetic pressure 
does not affect the hydrostatic equilibrium, except in the innermost regions, i.e. 
r <^ 1 kpc (see Gongalves & Friaga 1998 for a more detailed discussion). 
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